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^ . A simple, often invoked, regularization scheme of quantum mechanical path in- 

tegrals in curved space is mode regularization: one expands fields into a Fourier 
series, performs calculations with only the first M modes, and at the end takes the 
limit M —>■ oo. This simple scheme does not manifestly preserve reparametrization 
invariance of the target manifold: particular noncovariant terms of order h 2 must be 
added to the action in order to maintain general coordinate invariance. Regulariza- 
tion by time slicing requires a different set of terms of order h 2 which can be derived 
from Weyl ordering of the Hamiltonian. With these counterterms both schemes give 
the same answers to all orders of loops. As a check we perform the three-loop calcu- 
lation of the trace anomaly in four dimensions in both schemes. We also present a 
diagrammatic proof of Matthews' theorem that phase space and configuration space 
path integrals are equal. 
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I Introduction 



Quantum mechanical path integrals in curved target spaces have various applica- 
tions. Two decades ago they were used to quantize collective coordinates of solitons 
M, and one decade ago to calculate chiral |2| and trace [^, f|] anomalies in quantum 
field theories. In the former case, the path integrals were regularized by expanding 
the fluctuations about the solitons into normal modes, and cutting the sum over 
modes off at a maximum energy or at a maximum number of modes In the 
trivial vacuum, on the other hand, two methods have been widely used: the time 
discretization methocf] which has been pioneered by Dirac and Feynman, and the 
mode truncation method already described in |J. In the literature, both methods 
have been used on a par. 

In this article we point out the relation and the differences between the two 
methods. Both the time discretization method and the mode truncation method 
are particular regularization schemes which define the path integral. One expects 
therefore that the results for physical objects, for example the transition element 
T = (xf,tf\Xi ,U) , can be obtained from path integrals whose actions differ at most 
by finite local counterterms. In quantum gauge field theories, one can fix such am- 
biguities by requiring that Ward identities are satisfied. Similarly, for theories in 
curved target space which are general coordinate invariant one can require that the 
transition element be a "bi-scalar" (a scalar in x^ and xf), but this fixes the ambigu- 
ities only up to covariant terms, namely a term proportional to the scalar curvature 
R. Only an experiment, for example the measurement of the trace anomaly, can fix 
the coefficient of this term. One can also view the Hamiltonian operator H as an 
observable, and by writing T = (xf, , U) = (ij! |exp(— ^H)\xf) where f3 = tf — U 
in Euclidean space, one can fix the ambiguities in T by requiring that it satisfies 
the Schroedinger equation with H as Hamiltonian. If H itself contains no term 
proportional to R, then the action will contain a term — |i? as we will show, and 

T will have an R term in the exponent with coefficient —■h- If H contains a term 
aR, the coefficient of R in T will be — ( — + a). In general different regularization 
schemes can lead to different finite local counterterms of any order in h, but finite 
local counterterms of order h and h 2 only are generated in the path integral action if 
one considers different operator orderings in the Hamiltonian. A Hamiltonian which 
is a target space scalar^ is 

H = ^g~^(x)p^(x)g^(x)p u g-^(x) (1) 

but changing the order of the operators will in general destroy the invariance under 
target space diffeomorphisms. If the ambiguities due to using different regularization 

1 Time discretization in the sector with solitons is complicated because the canonical mo- 
menta (of the non-zero- modes) satisfy the equal-time commutation relation [tt(x, i), tp(y, t)] = 
— ih[6(x — y) — (f>' so i(x)(p' sol (y) /M so i] where 4>soi(x) is the classical soliton solution with mass M so i. 

2 The generator for target space diffeomorphisms x^ — ► x^ + is a sum of an orbital part 

G° rb (0 = sk + e(£)Pv) and a spin part = {28^9^ - ^A^afr Closure 

of the commutator algebra fixes G spm (£). An example of a target space diffeomorphism scalar 
is (|l|), namely [G orb (£) + G spm '(£) , H] = 0, as one may check by an explicit (tedious) calculation. 
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schemes correspond to the ambiguities due to different operator orderings, the terms 
of order h 3 and higher in the Hamiltonian should be the same in all cases. This 
conclusion is corroborated by the observation that in these non-linear sigma models, 
iV-loop graphs are convergent by power counting for N > 3. 

In the time discretization method, a very clear connection between H and the 
transition element exists if one uses phase space path integrals (Feynman's approach) 
instead of configuration space path integrals (Dirac's approach), and rewrites the 
Hamiltonian in Weyl ordered form. Then one may replace the Weyl ordered operator 
H-w(p,x) in the kernel by the corresponding function at the midpoint (Berezin's 
theorem) 

dp (x 1 \exp(--H w )\p)(p\x 2 ) -> exp(--H(p, 1 - 2 )j J dp (x 1 \p)(p\x 2 ). (2) 

The Weyl ordering of H in (1) leads to the following local finite correction to the 
naive Hamiltonian H = ^g^p^Pu to be used in the path integral 

AV w = T ^(R + g^T^T ua n, (3) 

where T^ u a are the Christoffel symbols.]] This counterterm has been extensively 
discussed and used in the literature 0. 

Both in the time discretization and in the mode truncation method, one writes 
the paths as x^(t) = x% (t) +q fM (r) for — 1 < r < 0, where x% (t) is a solution of the 
free field equations with the correct boundary conditions (x^ g (r) = x'i + t(x^ —xf)) 
and one expands the quantum fluctuations ? m (t) into a complete set of eigenf unctions 
of the free field equations with vanishing boundary values and integrates over a finite 
number of coefficients q^. The difference between the time discretization method 
and mode regularization is that in the latter one uses the continuum action and 
naive continuum Feynman rules for a finite number of modes and one chooses a 
measure for the integrals over these modes which is usually normalized such that 
it reproduces the standard result for the free particle if the interaction part of the 
action, Si nt , is set to zero. In the time discretization method, on the other hand, one 
uses a discretized action whose discretized Feynman rules are derived (in explicit 
form ||) from the Hamiltonian starting point. 

Up to this point we have made fairly obvious statements. Both regularization 
schemes are well-defined, and this is how they are used in the literature. However, 
it is already clear that this cannot be the whole truth, because with time slicing one 
has the nontrivial AVyy in (|3|), whereas with mode truncation such a AV^ seems at 
first to be absent. The precise way in which also mode truncation leads to a AV 
will be presented in section [TI[ We will see that one needs the following counterterm 
for mode regularization 

AV MR = j(r~ \g^g aP g,8T^T u ^ . (4) 

3 Our conventions for the Riemann and Ricci tensor are: R piiv a — dpT^ + T pT a Y pu T — (p <-> /i); 
R plJ = Rpuv 1 *- Hence, at the linearized level, — \ (d^d^h — dph v — d u hp + Uh^J) where 

h = rfgp U and = rf p d p hp V . 
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As a check we shall obtain the correct trace anomaly in four dimensions from a three- 
loop calculation. In Riemann normal coordinates, the Christoffel symbols V ^ v p in 
AV vanish at the origin of the coordinate system, but at the three- loop level one 
finds a contribution from AV by expanding each Christoffel symbol into q T d T T fll/ p 
and contracting the two quantum fields. Since <9 T r MJ , p contains a part proportional 
to the Riemann tensor, it is clear that AV yields a nonvanishing contribution. Hence 
dropping AV^ in Riemann normal coordinates is incorrect. 

If one eliminates the momenta by integrating over them in the time discretized 
path integral, one obtains N factors (detg Ml/ )~2 at the N midpoint coordinates 
!(<7fc+i + Ik)- Exponentiating these in the familiar Faddeev- Popov way one obtains 
what we have called "Lee- Yang ghosts" , namely commuting a p and anticommuting 
U L and & ghost fields ||. Phase space path integrals are free from ambiguities 
because they are finite. Clearly, upon transition from phase space to configuration 
space, ambiguities are created; technically this is due to the fact that the momenta 
are replaced by the q 's and ghosts, each of which introduces divergences and hence 
ambiguities into the theory. 

In the phase space approach, the vertices are different (for example \p^ h \x)p v 
instead of ^x^g^^x)^), as well as the propagators {{p^p y ) is not proportional to 
9fiagu/3(q a q t3 )), but the transition element should be the same (Matthews' theorem). 
In section |Hj we present a graphical proof. 

In section |V| we draw conclusions and we show that the coefficient of the R term 
in AV is scheme independent and equals |. 



II Mode regularization 



We now describe how one can define mode regularized path integrals in curved 
space. Ideally, one would like to derive mode regularization from first principles, i.e., 
starting from the transition amplitude defined as a matrix element of the evolution 
operator acting in the Hilbert space of physical states: (xy |exp(— j-H) \x% ). However, 
this derivation looks quite difficult, so that we prefer to take a more pragmatic course 
of action, and attempt a direct definition of the mode regularized path integral. This 
definition will be supplemented by certain consistency requirements which we specify 
later on. 

The transition amplitude can formally be written as follows 



x(0)=Xf r ^ 



(3% S . 



(5) 



{xp tf\x1 , ti) = I Vx exp 

•> x(-l) = x i 

S = J dr ^g^{x){x»x v + a^a v + 
Vx= [] d D x(T)d D a(T)d D b(r)d D c{T). 

-1<T<0 

We have shifted and rescaled the time parameter t — tf = /3t, and since all terms in 
the action only depend on j3h, we set % = 1 from now on. Note that (3 then counts 
the number of loops. We will evaluate the path integral in a perturbative expansion 
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in (3 and in the coordinate displacements £ M about the final point: £ M = xf — x'j. 
Formally the path integral is a scalar since the action is a scalar and the ghost fields 
make up a scalar measure on the space of paths.0 Mode regularization will destroy 
this formal covariance. However, we will see that covariance can be recovered by 
adding a suitable noncovariant counterterm AV to the action S. 
We start parametrizing 

^(T) = <(r)+^(r), (6) 

where x^ g {r) is a background trajectory and q^ij) the quantum fluctuations. The 
background trajectory is taken to satisfy the free field equations of motion and is 

A 1 
f 

enforcing the boundary conditions 



a function linear in r connecting xf to x*i in the chosen coordinate system, thus 



^(r)=x;-fr; ^ = x^-xff. (7) 

Then the quantum fields q^{r) should vanish at the time boundaries and can be 
expanded into sines. For the Lee- Yang ghosts we use the same Fourier expansion; 
this may be considered as part of our definition of mode regularization.0 Hence 

oo oo 

^( T ) = E <sin(7rmr); a^(r) = ^ a^sin(Trmr); 

m=l m=l 

oo oo 

= E Csm(vrmr); c^r) = £ Csin(vrmr). (8) 

m=l m=l 

At this point the formal measure T>x can be defined in terms of integration over 
the Fourier coefficients 

oo D 

Vx = VqVaVbVc = A ]J J] dq^da^db^, (9) 

m=l (1=1 

which fixes the path integral for a free particle up to the constant A 

J f)x exp —~aQ — A) Q = /_° dr ^M^T + o?a v + Wc v ). (10) 

The constant A will be fixed later on from a consistency requirement. Note that 
limiting the integration over the number of modes to a finite number M gives a nat- 
ural regularization of the path integral. This regularization resolves the ambiguities 
that show up in the continuum limit. 



4 The factors \\{det g{xi)) 1 ' 2 dx 1 ^ in the discretized path integral are target space scalars and 
exponentiating them leads to ghosts. 

5 Another argument to justify that the ghosts should be expanded into sines is that the classical 
solutions of their field equations are a M = = c M = 0, and that the quantum fluctuations do not 
modify the boundary conditions of the classical solutions. In J§| it was shown that the results for 
the transition amplitude do not change if one uses cosines for the ghosts. 
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Now we expand the action about the final point x^ and obtain 

S = 5*2 + Si n t = 5*2 + 5*3 + 5*4 + . . . 

where 

Ss = f dr l -d a g, v {q a ~ + ?<t + + 6V - 2q^ 



(12) 



S 4 = f dr -d a d p g^{q a qP + Cfr 2 - 2q a ^r)(^ u + q»q v + a»a v + &"c" - 2q^ u ). 
J — i 4 

All geometrical quantities, like g^ u or d a g^ u , are evaluated at the final point xf, but 
for notational simplicity we do not exhibit this dependence. S 2 is taken as the free 
part and defines the propagators while S int gives the vertices as usual. Therefore, 
the quantum perturbative expansion reads: 



(xf,t f \x?,ti) = / (Vx) exp 



P. 



Ae 2 " 



(e mt ) 



Ab-**'*'*' ((i - ^ 3 - l -s, + ^ 2 sl) + otfi) 



(13) 



Aiming at a two-loop computation, we have kept only those terms contributing up 
to 0(j3), taking into account that £ M ~ 0(y/]3), as follows from the exponential 
appearing in the last line of (|1"3|). The propagators that follow from S 2 are given by 



(14) 



(<f(r)g»> = -W(x f )A( T ,a) 
(a»(ry(*)) = W(x f )"A(T,a) 
(&"(t)c») = -2W(x f )"A(T,a) 

where A is regulated by the mode cut-off introduced below (|i~0|) : 



M 



A(r,(7)= £ 



m=l 



TC 2 m 2 



sin irrriT sin nma 



and has as continuum value 

A(r, tr) = r(<r + l)0(r - a) + tr(r + l)0(tr - r) 



(15) 



(16) 



(A dot on the left(right) of A(r,a) indicates differentiation with respect to r(a).) 
Using standard Wick contractions, we computed: 



<-/3> 



(17) 



(g^C^ + g a/3 ^C - 2g^ a i v i 



1 1 

"iS 12 



24 



24 



'2/3 



S%) — dag^dpgxp 



—{g aP g^g Xp - ^g ap g^g px - 6# Q VV P 

96 



+Ag ap g^g vX + 4# a VV P ) 

+i (VVW + % a V A - g aX g^)Ce + (2g aX g Pp - g aP g Xp )eC 



48 
p 9o 



(19) 



This completes the calculation of the transition amplitude in the two-loop approxi- 
mation using the mode regularized path integral. 

At this point we should make contact with other schemes and test the consistency 
of our rules. To do that, we use our transition amplitude to obtain the time evolution 
of an arbitrary wave function 



/ d D Xi^Jg~{^}(x^ t f \x p , U)*{xi, U). (20) 



We need the factor yg(xj) because the transition element is formally a biscalar as 
we explained before, but then ty(xf,tf) is a scalar and hence also ^(xj, which in 
turn implies that the measure must be a scalar as well. 

Since the transition amplitude ([13]) is given in terms of an expansion around 
the final point (xf,tf), we Taylor expand also the wave function \I/(xj,tj) and the 



measure J g(xi) in eq. (P0|) about this point 

VfaU) = *(x f ,t f )-pdMx f ,t f )+ed^(x f ,tf) + ±eCd»d^(x f ,t f ) + 0(pi: 



g ( Xi ) = ^g( Xf )^ + eT m a +-ec(d^ va a +T m a r v /)+o(e)). (21) 

We then perform the integration over d D x fl = d ^ in (|20| ) and match the various 
terms. The leading term fixes A 

* = A(2irf3)%^ -> A = (2ir(3)~%, (22) 

and the terms of order (3 give 



-d t + l - V 2 * + ±m - ^g pv g ap g^d,g ai d u g^ + ^g^g^d.g^g^ 



0. 



(23) 

This last equation means that the wave function \1/ satisfies the following Schroe- 
dinger equation at the final point (x p ,tf): 

- d t V = (H + AV e// )tf = -^V 2 ^ + AV e ffV, (24) 



where the effective potential is given by 



= ~R - ^IV/IV* + ^g^g^d^d^s (25) 

= -\R+^ v g aP g,sT^Vj. 

Clearly, to obtain the "free" Hamiltonian H from a path integral, one should sub- 
tract the potential AV e ff, and thus use in (||) the following classical action 

-i£ = /V (-± g ^ x )x»x» + f3AV eff (x)y (26) 

With this counterterm mode regularization gives the same results as the time dis- 



cretization scheme described in the introduction. Note that (|25| ) is different from 
the counterterm for time discretization. The difference is the last term in the second 
line of 

Finally, as a non-trivial test, which was one of our motivations to carry out the 
present research, we have computed the trace anomaly for a real conformal scalar 
field in four dimensions using mode regularization and the newly found counterterm. 
Anomalies in Fujikawa's approach are given by the regulated trace of the Jacobian 
J of the symmetry transformation. 

An = limTr(Je-t^). (27) 

A theory of consistent regulators TZ exists and it gives us for conformal scalar fields 
the Hamiltonian (|l|) minus an improvement potential % %n^\\ R- The Jacobian for 
trace anomalies is proportional to unity. In mode regularization ( ^7|) translates into 



the evaluation of the expectation value of the Jacobian with respect to the path 
integral based on the formal transition element (|5]) including the just calculated 
counterterms and the improvement term. Taking the trace of (|5]) means first putting 



Xf = Xi which gives that x^g = Xf and then integrating over Xf with d D Xf^Jg(xf). 

As the measure of (||) includes the constant A = (27r/?)~"? ; the /^-independent term 
of the r.h.s. of ([27]) is given by the (D/2 + l)-loop term in the path integral. 
Expanding the metric in Riemann normal coordinates we find the first nonvanishing 
contributions at the two-loop level: graphs with the topology of the number 8 and 
the background values of the various potential terms. From (|18|) with £ M = one 
obtains^ 

+ . ^ ( - w if4W( W i - (g -; )/(g - i) v 



6 V 4 

The result is the same as obtained from time slicing and yields the correct trace 
anomaly in D = 2 dimensions. 
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In equations (p8|), (p9|), (30), ( |3l| ) only the topology of the graphs is indicated, e.g., the figure 
in stands for all possible graphs of that shape, including ghost loops. 
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At the three-loop level we find three regular graphs 




72 
1 

72 



-m 2 
i 



-Is- 2 

6 ^ J ' 
-I/? 2 



(29) 
(30) 



as well clS cL graph coming from the potential terms 




-\(BKf 



1-(D-2)/(D-1) v2r _ Lr2 



16 



72 



(32) 



There are only two differences with respect to time slicing: with time slicing, the 
factor — | in (j30|) becomes — |, and the factor — ^ in (|32|) becomes — -h- These 



two modifications lead to the same final expression. Adding all contributions of 
connected and disconnected graphs we find the correct result 



An(Weyl)(spinO,D = 4) = J — ^ Jgjxj a (x) — (i?^ Ap - R% - V 2 i?) . (33) 

where <j(x) is the arbitrary function appearing in the Jacobian: J = a(x)5 D (x — y). 
This shows that mode regularization works after all.[] The fact that the right answer 
is obtained suggests that no new counterterms are needed in this scheme. A complete 
three-loop calculation in arbitrary coordinates and with a non-vanishing £ M could 
be used to test eq. (|(J) at order f3 2 . As we noted in the introduction, however, all 
three- and higher- loop graphs are power-counting convergent. This means that at 
these orders in h any consistent regularization scheme will yield the same answer. 
We therefore claim that with the presently found counterterm, the mode regularized 
path integral is consistent to all orders. This suggests that the difference is indeed 
due to operator-orderings. It would be interesting to justify mode regularization 
from first principles. 

Before closing this section, we comment on how loop integrations in Feynman 
diagrams are done. Mode cut-off allows one to disentangle ambiguities that appear 
in the continuum limit of certain integrals over the A's. Resorting to the mode 
regulated expression for A, one can use partial integration and take boundary terms 
into account if they are non- vanishing. Using partial integration repeatedly, one gets 
expressions containing only A's, A"s and *A's which are unambiguous in the con- 
tinuum limit, and computes them there. Useful identities obtained at the regulated 

7 In Q the noncovariant part of the counterterm (|2^) was missed. In Riemann normal coordi- 
nates this term affects only the coefficient of the R 2 ^ v xp term in the four dimensional trace anomaly, 
which was erroneously calculated. 
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level from (|T5| ) are 



'A'(r,T) + "A(T,T)=d T ('A(r,r)), 

"A(r,a) = A"(r,a). (34) 

The first identity was used to compute (|28|). Here is a list of integrals that are easily 
computed using partial integration?] and whose values are different from the time 
discretization method: 

h = Jdrr ('A' + " A) | CT=T = 0, 

h = JJ drda ('A) (A«) ('A') = -±, (35) 
Js = // rfrrfa A («A) (A*) ('A') = ^. 

Time discretization would give Ji = — |, J 2 = — | and -^3 = 36o'0 

The non-vanishing value of l\ in time discretization is needed to compensate the 
explicit factors of g^ l ^(xi) appearing in the measure, see (|36|) , and together with 
I2 it is responsible for the different counterterms required in the two regularization 
schemes. Finally, -Z3 leads to the different values for the coefficient of R 2 ^ u \ p in 
eq. ( |30|) in the two schemes. The counterterms with two Christoffel symbols are also 
different in both schemes but the final result for the transition element (and hence 
the trace anomaly) is the same. This confirms our approach to mode regularization. 



Ill Phase space path integrals in curved space 
and Matthews' theorem 

In the phase space approach to path integrals in curved space the vertices and 
propagators are different from those used in the configuration space approach. For 
example, the leading terms in the action are \p lJ ,g^ v {x)p v and ^x^g^x)^ , respec- 
tively, and the (g M (r)g^ '(er)) propagator contains a Dirac delta function but the 
(p^(r)p u (a)) propagator is completely regular. Although the vertices and propa- 
gators differ, the result for the transition element should be the same (Matthews' 
theorem^) [|K| , and the way this comes about is due to a new kind of ghosts |3|, £| . 

8 To obtain J 2 , use (A , )(*A*) = ^d T (A*) 2 , partially integrate, use "A = A", and then use 
(A")(A*) 2 = ±<9 CT (A') 3 . To obtain I 3 , write the integrand as A(*A)i<9 r (A') 2 , partially integrate, 
and use then that A(**A)(A*) 2 = A(A")(A*) 2 = A|d T (A*) 3 . Then partially integrate once 
more. 

9 In the time discretization scheme 5(a — r) is a Kronecker delta function and 9(0) = \. This 
leads to a fully consistent scheme as shown in M . For example, J J da dr 5 (a — r)0(a — t)6(<j—t) = 
J J da dr 5(a — r)0(a — t)6(t — a) = j whereas mode regularization would give | and \ ~ \ = \ 
respectively. Using *A(r, a) = a + 9(t — a), "A'(r, a) = 1 — S(t — a) and "A = S(t — a), the 
results below (|35| ) follow. 

10 Matthews' theorem holds for connected graphs and originally only applied to nonderivative 
interactions. 
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The difference proportional to 8(t — cr) between the pp and qq propagators turns out 
to be the propagator of those ghosts which we have called Lee- Yang ghosts after Lee 
and Yang, who where the first to point out that in the deformed harmonic oscillator 
(using arbitrary coordinates in one dimension) one needs extra 6(0) terms in the 
action or Hamiltonian |jTIfl . The modern way to deal with these ill-defined objects 
5(0) is to introduce new ghosts, just as in the familiar case of Faddeev-Popov ghosts 
for gauge theories. From the formal definition of the path integral, it is clear that 
Matthews' theorem should hold but as always it is desirable to check this by 
a nontrivial calculation. Moreover, such a calculation in turn may reveal a relation 
between the Feynman graphs in both approaches which can lead to a diagrammatic 
proof of Matthews' theorem. With this in mind we calculate the two-loop correction 
to the transition element (xj, 0\xf , —(3), both using the phase space approach and 
the configuration space approach to path integrals. 

Using time discretization (see footnote |9|), the transition element is given in both 
phase- and configuration-space by 



(x f \exp(~k)\xi) 



9{xj) 



1/4 



(2n(3h) D / 



(36) 



D l 2 is the usual Feynman measure and (xf\ and \xi) are both normal- 

1/4 



where (2ir/3h) 

ized as J d D x^J g(x)\x) (x\ = 1. The factor 
from this normalization.^ Furthermore, r 



9(*f) 



X 



bg 



V) 



_ t_ 

P 



in the measure is a direct result 
and x^(r) = x^ g (r) + q^(r), where 



,1C ■ 



is the background solution of the free equations which 



satisfies the correct boundary conditions, and q^(r) are the quantum fluctuations. 
The phase space interactions and propagators are given by 



1 npha 



1 

"j3h 
-ph 



dr - 
-i 2 



grid 



o 



dr AVi 



w , 



(37) 



(q^r)p u (a)) = -z/3^[r + 9(a - r)] ~ zA*, (38) 
(qW(a)) = -Phsr{xf) A(r, a) ~ A, 

where we recall 

A(r, a) = r{a + l)0(r - a) + a(r + 1)9 (cr - r). (39) 

We have rescaled p M (cr) into 4p M ((r) in order that all propagators have a factor (3% 
and all vertices a factor (/3^.) _1 . The (3 and h appear always together and count 

11 Inserting iV sets of p-eigenstates and N—l sets of x-eigenstates gives a factor [g(x f)g(xi)} 1 ^ , 
but the phase space path integration over the p's and the (?'s gives an extra factor g(xf)i because 
there is one more p than g's. This extra factor is taken at the point Xf because we expanded 
the metric around this point. In the configuration space path integral, one integrates over the 
p's, exponentiates the factors g(x) by introducing ghosts, and then produces factors of g(x) by 
integrating over the ghosts, with the same final result. 
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the number of loops. The configuration space interactions and propagators are as 
in section [II], though regulated as in footnote |9| 

- \S^ S = ~ £ dr ~ (g, v {x) - g^Xf)) (x»x» + + b»c v ) - (3h £ dr AV W ; 



(<f(r)g> 
(9"(t)«T((7 
(<f(r)?> 
(a^(r)a u (a 
(^(r)cV 



-(3g^(x f )(r + e(a-T))^ 
-2(3g^(x f )5(r - a) ~ 2"A 



A* 
•A* 



(40) 



From *A* = 1 — 5{r — a) and "A = S(r — a) it is clear that the sum of the (qq) 
and the ghost propagators is equal to minus the (pp) propagator. 



+ • • 



. «/VA/\A/WW» 



(41) 



Further, the (qp) and the (qq) propagators are equal up to a factor i (note that 
p = ix is a field equation). 



x i 



(42) 



We are now ready to evaluate two-loop diagrams. The result for the £ M 
independent part is given in figure |T] In the phase space approach, one finds deriva- 
tives of the inverse metric g^ v which we have converted to derivatives of the metric 
g^y for purpose of comparison. The relations A, B, F, and G follow immediately 
from (f|l]) and (]42). Considering the identities C, D and E one must note that the 
phase space expression comes from both one- and two-vertex graphs. The identities 
then follow from rewriting fl4lD as 



«A/WV\A/W» 



and replacing a "ghost" propagator on the phase space side by its value, (—f3K) 
times a delta-function, which pinches the two vertices. It is clear that the two- 
loop corrections in configuration space agree with those in phase space. We have 
also checked the ^ dependent terms in the transition element and found complete 
agreement with ([13|). As, once again, all configuration space graphs at the three- 
and higher loop level are unambiguous, the two-loop calculation suffices to show the 
equivalence between the phase- and configuration-space path integrals. 



IV Conclusions 



The two regularization schemes considered in this article, mode regularization 
and time slicing, give the same answers at all orders of loops provided one adds 
specific order h 2 counterterms to the action in the configuration space path inte- 
gral. These terms are proportional to the curvature R but also to products of two 



12 



Phase space 



Configuration space 




Figure 1.: Two-loop diagrammatic identities using time slicing. Each corresponds 
to the coefficient of a particular metric structure present in the two-loop parti- 
tion function; A: -\g^n 9liV , B: -\g^g^d a d p g pu , C: -\g a ? g» p g va d a g^d p g pa , 
D , -±g^gPpg™d a g^dpg pa , E: -\ g af> g^g u °d a g^d p g pa , F: -\g^ ! <T 'g pc d a g^g i 
andG: -\g^g fiv g^d^d^. 

Christoffel symbols and the latter are different for the two schemes. These extra 
terms of order h 2 follow in the time slicing method from rewriting the Hamiltonian 
in Weyl-ordered form, but for mode regularization we could only determine them by 
requiring that the Schroedinger equation be satisfied. (This way of fixing the extra 
terms can also be applied to time slicing). As a nontrivial check that these extra 
terms are indeed needed (and also a check on their explicit form) we evaluated the 
three-loop contributions to the trace anomaly in both schemes and indeed found the 
same (correct) answer. 

We also considered phase space path integrals. Their loop graphs contain no 
divergences. We gave a diagrammatic proof of Matthews' theorem that phase and 
configuration space path integrals give the same transition element. A key obser- 
vation was that the difference between the phase space propagator (p(a)p(r)) and 
the configuration space propagator (q(o')q(T)) is equal to the propagator for the 

Lee- Yang ghosts. The latter result when one exponentiates the factors \J g(x) which 
are produced by integrating out the momenta. 

It remains to explain why the coefficient of the R term in AV is equal to | in all 
"reasonable" regularization schemes. Let us call a scheme reasonable if it satisfies 
Matthews' theorem. Then the first two graphs in figure 1 should give the same 
result in configuration space as in phase space. The equality for graph B is always 
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satisfied, see ( p8|) and (^3), but for graph A one finds the condition 



dr A( 



r, r 



A( 



r, r 



*A*| 



r, a 



"A( 



(43) 



Both time slicing and mode regularization satisfy this condition. In all such rea- 

g fMU (x)x fl x u will 

then contain a two-loop contribution just as the phase space calculation would 
give. The complete phase space Hamiltonian fl3"T|) contains, however, in addition the 
counterterm (| 

any reasonable regularization scheme and the R term in the action to be used in the 
configuration space path integral has coefficient 



sonable schemes the transition element with the naive action S — \ 

h 2 



The total R term in the Hamiltonian element is then — j^R for 
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